IN THE SPECIFICATION; 



(1) Kindly replace the paragr^h beginning on page 10, line 4 with the following 



paragraph: 



—The geometry of the nets CI, C2, C3 is described hierarchically and captured to a level of 
detaiTiTeeded-tOLd eiermine an accurate solution. For example, in FIGURE 1 , the net CI and the net 
C2 are captured with relatively detailed geometric descriptions. However, the net 
captured with relatively coarse geometric descriptions.— 





(2) . Kindly repl^^^^ttj^^aragraph beginning on page 7, line 8 with the following 
paragraphs: 

—FIGURE 2A illustrates a block diagram of a system for determining a capacitance of an 
integrated circuit constructed according to the principles of the present invention; 

FIGURE 2B illustrates a flow chart of a method of determining a charge distribution for a 
net constructed according to the principles of the present invention;— 



(3) Kindly replace the^p^:^raph beginning on page 11, line 5 with the following 
paragraph: 

— Tuming now to FIGURE 2A, illustrated is a block diagram of a system for determining a 
capacitance of an IC, generally designated 200, constructed according to the principles of the present 
invention. The system 200 comprises an integral equation formulator 210, a charge variation 
function generator 220, and a conductive geometry generator 230. The integral equation formulator 
210 determines the total capacitance of all the nets CI, C2, C3 within the IC— 
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(4) - ^ Kindly replace Aerparagraphs beginning on page 1 3 , line^20 and ending on page 24, 



)^agraphs beginning on pa^i 
^e 10 with the following paragraphs: 



"In one embodiment of the present invention, the characteristic geometry function for the 
set of rief surfaces; Xr*, is a predefined function. In an alternate embodiment of the present invention, 
Xr is calculated from information related to the. capacitive structure of the IC. Xr is scalar- valued 
function that has a n-dimensional position as thelnput. To illustrate the form of Xr, suppose the 

A 

geometry to be described is the unit square in the jcy-plane! The fprm of for this geometry then 
would be: - 



X^~s{x)s{y)5{z\ (1) 



% and^( 



where s(x)=l for 0<x< 1/and s(x) is zero elsewhere, and d(z) is the standard Dirac delta function. 
Those skilled in the art are familiar with the properties arid use of Dirac delta functions. 

The function f, which represents the charge distribution v^thout regard to geometry, is 
evaluated over R, the set of surfaces of the nets CI, C2, C3. The surface integrals over R are 
reformulated as volume integrals via the relation: 



i/=lA» (2) 



where V is the volume that the net comprises. 

To solve the iterative linear solution, the integral equation formulator 210 uses an initial 
guess g, for the charge distribution, and a starting subdivision of the geometry of the net. The 
integral equation formulator 210 also creates the projection matrix Pr from Xr, 



m 



(3) 



The bracketed expression above is the standard inner product for functions: 
^ {f,g)=\fg (4) 

Also in this equation, P,^ is a primitive projection matrix where the (ij) entry of/^ = ^IJjl/^ . 
Finally the IjS represent the choice of orthogonal polynomials. In one embodiment of the present 
invention, Legendre polynomials form the system of orthogonal polynomials. One skilled in the 
pertinent art is familiar with and the use of Legendre polynomials. Also, P^'r will represent that part 
of Pr that contains information oftiie geometry of subdivision b. 

Using the, initial guess g for the^charge-distribution and the starting subdivision of the 
geometry of the net, the integral equation formulator>2 1 0 determines charge distributions for a given 

potential on the nets C 1 , C2, C3 . From the charge distributions, the integral equation formulator 2 1 0 

'\ 

determines the capacitance of the n^ts 1 1 , C2, C3. 

In the determination of the charge distributions,'' the iterative linear solution uses the 
following fundamental equation that relates the charge and the potential: 



i^{r)=lG{ry)p(r')dr' ' (5) 



where G (r, r') is the Green's function for the geometry of the IC and" p(r'Xis the charge, density. 
Those skilled in the art are familiar with the properties and use of Green's functions. 



• 



In one embodiment of the present invention, the fundamental equation relating the charge 
and the potential is expressed as: 

yr = MP J (6) 

where the v|/ is a potential distribution and f represents the charge distribution without regard to the 

geometry. From this equation, the integral equation formulator 210 can determine the charge 

\ 

distribution function f. 

In one embodiment of the present invention, the integral equation formulator 210 uses a 
Krylov method based on a variation of a Generalized Minimal Residual Method (GMRES) for the 
iterative linear solution. Backgrouridnnformation conceming GMRES is discussed in GMRES: A 
Generalized Minimal Residual Algorithin for Solving Nonsvmmetric Linear Svstems , by Y. Saad 
and M. H. Shultz, SIAM Journal on .Seieritific"and Statistical Computing, 7(3):856-69, 1986 
(incorporated herein by reference). \ 

GMRES is^used to solve for the'charge distribution function f in the following equation: 

l^/?/-{^olL=^ 

or, equivalently: 

P,MPJ=P,ilf,, (8) •• . ■ 

where Mis an operator that converts a charge density to a potential distribution and \|/o is a potential 
which is 1 over an individual net and 0 elsewhere. Those skilled in the pertinent art are familiar with 
the application of GMRES. Nevertheless,, below is a description of the, process both to give a 



complete picture of the iterative linear solution process and to show the advantages of decoupling 
the cljferge variation from the geometry. 

The integral equation formulato^2 10 starts iterative linear solution with an initial guess g for 
the charge distribution and an initial geometry. The integral equation formulator 2 1 0 also determines 
X^'r, which represents the geometry of the box b in the subdivision, and ?\, which is the part of the 
projection matrix which covers the box b (see FIGURE 3 A for more information concerning 
geometry subdivisions and boxes). From this information, the integral equation formulator 210 
computes a potential y based upon the initial charge distribution guess g and the current subdivision. 

The integral equation formulator 210 also determines the desired potential \\fo:^ 

P,Wo-PsO (9) 

where o is the constant function 1 oyer anhidividu&l^net and Pg is defined as Pr over an individual 
net S. The integral equation formulator 210 then^determines the difference between the desired 
potential \|/o and the potential y. This difference, v|/o-\|/, becomes first charge variation function f,. 
Note that fj is not a representation of the charge distribution described by f fj is a fimction used to 
modify the charge distribution function f to reach a desired resolution. 

s 

The integral equation formulator 210 determines the ratio y/p, where y = llfJlR and p = lk|/oltR. 
If the ratio is sufficiently small, the iterative linear solution is complete. 

If the ratio is unacceptably large, the integral equation formulator 210 normalizes f, and 
proceeds with the linear iterative solution. The integral equation formulator 210 then invokes the 
charge variation function generator 220 to create a charge variation function which refines the 
description of the charge distribution. At the beginning of the iterative process the charge variation 
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function generator 220 creates a second charge variation function In general, the formula for the 
creation of the j+l* charge variation function from the j'^ charge variation function is: 

/y., = MP'jj . (10) 

Where M is a^ operator that converts a charge density to a potential distribution. The (ij) entry of 
the matrix representation of M is: ^ ^ 

M''= \^[l^rjG{ry)dr'dr (11) 

Where 1**^ is the i'*" moment in box b. Note that the operation of M involves all boxes, not just box 
b. G(r,r') is the Green's function foFthe set of net surfaces R. Those skilled in the art are familiar 
with the properties and use of Green's functions. 

The charge variation function "generator performs the operation M using a variant of the Fast 
Multipole Method (FMM) algorithm called the Fast Distributed Method (FDM). The FDM differs 
from the FMM in several areas. First, the FDM omits the FMM's initial step of computing multipole 
representations from point charges. FDM omits this step because the input is already in terms of a 
charge distribution. Second, the FDM omits FMM's direct point-to-point interactions. Instead, the 
FDM uses a source charge distribution in a cube. The potential distributions are calculated in all 
neighboring cubes, including the source cube itself. 

Third, the FDM omits FMM's final step of interpolating point potentials from the local 
expansions. The FDM omits the last step, because the desired output is a potential distribution. 
Finally, FDM uses Legendre polynomial expansions for both charge distributions and local 
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• 



expansions. The FMM uses multipole expansions for representing charge distributions and local 
\ expansions for representing potential distributions. 

The FMM algorithm, Legendre polynomial expansions and their use are well known to those 
skilled in the pertinent art. Background information concerning numerical analysis and capacitance 
calculations is discussed in Introduction to Numerical Analysis , by J. Soer and R. Bulirsh, Springer- 
Verlag 1 979 and in Preconditioned, Adaptive, Multipole- Accelerated Iterative Methods for Three- 
Dimensional First-Kind Integral Equations of Potential Theory , by K. Nabors, et al , SI AM Joumal 
on Scientific and Statistical Computing, 15(3):7 13-735, May 1994 (both incorporated herein by 
reference). 

The use of the FDM results in a significant savings in time over both the analytical 
calculation using the explicit form of M^'and the FMM. The savings are due to the omission of 
several FMM steps and the omission of FMM' s direct point-to-point interactions. 

Once the charge variation function generator 220 has created the charge variation function 
fj+i, the integral equation formulator 210 uses this new function to create the (i, j) entry of a 
Hessenberg matrix: 

^ = (12) 

foralli<j. 

The integral equation formulator 210 uses the Hessenberg matrix to solve for the new coefficients 
of the charge variation functions to be used to further refine the charge distribution fimction f 
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After the entries in the Hessenberg matrix are computed, the integral equation formulator 210 
then orthogonahzes fj+j from the other fjS. The charge variation fiinction generator 220 then creates 
the Ctlj) entry of the Hessenberg matrix by the expression: 



This expression is the standard norm of functions defined by: 



After the creation of this last entry of the hessenberg matrix, the integral equation formulator 210 
then normalizes fj+j. The orthogonalizatioh and the normalization of fj+j create orthonormal charge 
variation ftinctions. As a result, the integral equation formulator 210 is able to obtam the correct 
charge distribution function f very quickly. 

The integral equation formulator 210 then solves' for x, the coefficient vector for the charge 
variations, using least-squares on the equation: \ 

Hx^ye^ (15) 

where H is the Hessenberg matrix generated by the integral equation formulator 210 as explained 
above. is the unit vector with the form (1,0,0,0,0, . . .)(with j number of zeros). The vector x 
is the coefficient vector for the charge distributions and y is the norm of the first charge variation 
fiinction. 7 is defined as: 

r = ||/.IL- 06) 
-9- 




f ' Since the charge variation and the geometry are decoupled, the representation of the will not 
become very large as compared to the previous methods. Using the residual of the least squares, the 
integral equation formulator 210 computes the factor r as the 2-norm of the least squares residual. 

The integral equation formulator 210 then generates a new guess for the charge distribution 
function f using the equation: 

/=i 

The integral equation formulator 210 then determines the ratio r/p. Where r is the 2-norm 
of the least squares residual r and p is^he norm of the^lesired potential (defined earlier as p = Ikj/ollR). 
If the ratio is within acceptable limits, the charge distribution has converged and the integral equation 



formulator 210 termmates the iterative linear solution. 

\ 




If the ratio is not within acceptable limits, the integral equation formulator 210 uses the new 
guess for the charge distribution function f, calculated above, as the starting point for the next 
iteration of the iterative linear solution. However, before the integral equation formulator 210 
performs the next iteration, the integral equation formulator 210 must determine if the current 
subdivision requires a more refined approximation. 

The integral equation formulator 210 determines if the current subdivision b requires 
refinement by calculating a charge-geometry error e. In one embodiment of the^present invention, 
the integral equation formulator calculates the charge-geometry error e associated with box b in the 
subdivision using the following equation: 
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(18) 



i=k 



Where 1^ is the i**" moment basis function, f represents the distribution of the charge without any 
reference to the geometry of the IC and is the component of f along the i'^ moment for the charge 
function contained in the subdiyisiDiTb. The Simunation of Ij between k and N-1 represent the high 
order polynomials. Finally, N is the number of moments used in the expansion of f 

The integral equation formulator 2 1 0 then determines if the charge-geometry error is within 



acceptable limits using a charge-georiietry errpr criterion. The charge-geometry error criterion is 
defined as: 



R 



< 



(19) 



Where S is the minimal subdivision containing the selected net. £ is a small constant that determines 
the allowable tolerance of the criterion. Also, Hc^'Hr is the magnitude of the charge-geometry error 
in subdivision b. Note that the error threshold is relative to the charge distribution f only on the 
subdivision of the selected net.~ 




(5) Kindly insert the foll owing paragraphs after the paragraph that ends on 24. line 10: 
TuminglioWtafTOTJRF ?B,mth continued reference to FIGURE 2A, illustrated is a flow 
diagram of an embodiment of a method of determining a charge distHBlitteftJera, net, generally 
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designated 250, constructed according to the principles of the present invention. The method starts 
in a st^255 with an intent to determine a charge distribution. 

Arnriitial charge distribution and geometry are provided in a step 260. The initial charge 
distribution anH the initial geometry may be guesses and are used to start an iterative Hnear solution. 
The initial guess ibr the charge distribution may be designated g and the initial guess for the 
geometry may be a subdivision of the geometry of the net. 

After an initial\harge distribution and geometry are provided, a first charge variation 
function is then determined ik a step 265. The first charge variation function may be the difference 
between \\/q and \)/. In one embodiment, the first charge variation function may be determined by 
solving for v|/o and \|/ using Equations 6 and 9. 

After determining the first charge variation function, a determination is made if the charge 
distribution function is within an acceptable^mit in a first decisional step 270. In one embodiment, 
the acceptance of the charge distribution function may be within an acceptable limit if the ratio y/p 
is sufficiently small. If the charge distribution functris^n is within an acceptable limit, the method 250 
ends in a step 295. 

If the charge distribution function is not within anWceptable limit, the charge variation 
function is refined through an iterative linear method in a step 280. In one embodiment, f, may be 
normalized before proceeding with the linear iterative method. The^ethod 250 creates a charge 
variation function which refines the description of the charge distributioir^mploy ing Equations 1 0- 
17. 

After refining the charge variation function, a determination is made if the \Jiarge distribution 
function is within an acceptable limit in a step 282. In one embodiment, the chaige distribution 



mction is within an acceptable limit based on the ration r/p. If the ratio is within acceptable limits, 
the cn^ge distribution has converged and the method 250 continues to step 295 and ends. 

If th^^^etim-ge distribution function is not within acceptable limits, a determination is made 
if the geometry needs rfefmement in a third decisional step 287. In one embodiment, the method 250 
may employ Equations 1 8 anoH^ to determine if the geometry needs refinement. If the geometry 
does not need refinement, the method\;ontinues to step 280. If the geometry does need refinement, 
the geometry is subdivided into subdivision^ in a step 290. After subdividing the geometry, the 
method 250 continues to the step 280. 

While the methods disclosed herein have been'^scribed and shown with reference to 
particular steps performed in a particular order, it will be und^stood that these steps may be 
combined, subdivided or reordered to form an equivalent method wW^ut departing from the 
teachings of the present invention. Accordingly, unless specifically indicated herek\^the order and/or 
the grouping of the steps are not limitations of the present invention.-- 
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